Parallel Implementation of a Lattice Boltzmann Algorithm for the Electrostatic Plasma Turbulence 1 Lattice Boltzmann Algorithm for the Electrostatic Plasma Turbulence

A parallel version of a Lattice Boltzmann Equation algorithm , which simulates the electrostatic plasma turbulence, has been developed using the High Performance Fortran language. The algorithm evolves a system of particle populations on a discrete lattice and data-parallel implementation has been performed by a regular domain decomposition. System evolution requires both completely local and non-local operations, involving communication between processors. Communication phase has been minimized using local HPF procedures. EEciency tests of the parallel code, performed on a 9076 IBM-SP2 parallel computer , have given satisfactory results. A Lattice Boltzmann algorithm was developed to simulate the electrostatic turbulence in a magnetized plasma, within the framework of the two-uid model 1, 2]. The ion continuity, momentum and pressure equations are reproduced, with the electrons being described by the adiabatic response 3,4]. The Lattice Boltzmann Equation (LBE) method is a particularly promising numerical method to perform uid-dynamics high-resolution simulations 5]. The macroscopic dynamics is simulated starting from a system of particle populations neighboring sites, the coordinate index and D the lattice dimension. Particles which arrive at the same site undergo a collision, with the collision operator chosen in a way such to locally conserve the particle number and the total momentum. After the collision phase, each particle propagates toward one of the neighboring sites along the direction of its velocity. The LBE describes the evolution of such a particle system, assuming that, under the eeect of collisions, the system relaxes toward a local-equilibrium distribution function N eq j (r; t), and linearizing the collision operator around such distribution function. ?